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ABSTRACT 


In this paper a non-linear analysis of a M PLL by using the method 
of "Harmonic Balance" is presented. The particular M PLL considered has 
a low-pass filter and a band-pass filter in parallel . An analytic expres- 
sion for the relationship between .the input signal phase deviation and 
the phase error is determined for sinusoidal FM in the absence of noise. 

The expression is used to determine bounds on the proper operating region 
for the M PLL and to investigate the "Jump Phenomenon" previously observed.. 
From these results the proper modulation index, modulating frequency, etc. 
used for the design of a M PLL is determined. Data for the loop unlock 
boundary obtained from the theoretical expression are compared to data, 
obtained from analog computer simulations of the M PLL. 
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Ic INTRODUCTION 

The phase-lock loop (PLL) has become more widely used as an. FM , 
demodulator in the past several years. This is due to the superior 
noise rejection properties and threshold extension exhibited by the PLL. 
These properties are particularly useful in aerospace frequency modula- 
tion (FM) communication systems where receivers must operate in low 
signal-to-noise-ratio environments. However, the common low-pass PLL 
will, yield little or no threshold improvement over the conventional FM 
discriminator when the spectrum of the demodulated; FM signal contains in- 
formation at or near discreet frequencies in addition to the baseband 
signal. This is a characteristic of the typical IRIG spectrum. The 
Apollo unified S-band communication system is also such a system. . It has 
a spectrum (see Fig. 1) which consists of wide-band information at base- 
band and narrowband information on a subcarrier. In this case the low-pass 
PLL allows noise in the unused portion of the bandwidth tp enter the loop 
and thus lowers the loop signal-to-noise ratio. This problem is further 
compounded by poorer tracking performance of the low-pass, PLL when rela- 
tively high-energy subcarriers occur in the upper part of the passband. 

A solution to this problem .involves the use of additional filters in 
the loop in parallel with the low-pass filter [4]. This type of PLL is 
called a multifilter PLL (M PLL). .. The frequency response of such a M PLL 
matched to the above mentioned spectrum is shown in Fig.. 1 along with the 
conventional low-pass PLL. 

Recently there has been an .interest in a comparison between the 
multifilter phase feedback loop (M PFL), as used in precision ranging 
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systems [5], and the M PLL. It -has been shown [6] that the , parameters and 
filter responses of the M PFL and the M PLL are related. By using these 
relations the M PFL performance characteristics can be determined by the 
performance characteristics of .an equivalent M PLL. Therefore . the result 
of this .paper, can, be applied to the M PFL as well as to the M PLL. 

In this paper a. non-linear analysis of a M PLL is presented. The 
particular M PLL considered has a low-pass , filter and a band-pass filter 
in .parallel. An analytic expression for the relationship between the in- 
put signal phase deviation and the phase error is determined for sinusoidal 
FM in the absence of noise. The expression is used, to determine bounds on 
the proper operating region for the M PLL and to investigate the "Jump 
Phenomenon" previously observed [2] . From these results the proper modula- 
tion index,, modulating frequency, etc. used for the design of a M PLL is 
determined. Data for the loop unlock boundary obtained from the theoretical 
expression are compared to data obtained from analog computer simulations 
of the M PLL. 
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II. THE M PLL DIFFERENTIAL .EQUATION 


The Carrier Model of a typical phase-lock loop (PLL) is illustrated 
in Fig. 2. The phase detector generates an error voltage which is a func- 
tion of the phase difference (error) between the input signal and the VC0 
(Voltage Controlled Oscillator) output. This error voltageis then used 
(after filtering) to modify the VCO output so as to reduce the phase error. 
Thus the VCO attempts to "track" the input signal. Since the. VCO input 
is proportional to the VCO output frequency, the loop filter output is 
proportional to the input signal frequency (when the PLL is "tracking" 
the input signal). 

Viterbi [1] has reduced the Carrier Model to a baseband or phase model 
for the case in which the phase detector is simply an ideal multiplier 
(assuming, harmonics of the carrier frequency are filtered out). This 
model is shown in Fig. 3(a) where the. loop gain is assumed to be included 
in the loop filter. 

Assume that F(s) is a linear filter with a rational” transfer function 


given by 


F (s) 


, , 2 , m 

a + a,s + a„s + ....a s 
o 1 2 . m 


b + b,s + b„s^ + . . . .b s n 
o 1 2 n 


(la) 


where m <_ n. It can be shown [1] that the PLL operation is described by 
the differential equation. 

(b + b p + b p n )[|(t) - 8. (t) ] 

o I n i 

= -A (a + a,p + s ?.y.+ a p m ) ' s in <J> ( t ) 
o l m 


(lb) 
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where p is an. operator indicating differentiation. 

The phase model of a multi-filter phase-lock loop (M PLL) with a 
low-pass filter (LPF) and a band-pass filter (BPF) in parallel is shown 
in ; Fig. 3(b). This is just an extension of the model in Fig. 3(a). 

The LPF is necessary in- order that the M PLL track the carrier. The 
filter used in this study is one with a single pole transfer function given 
by 


F. (s) = 1 + - = 

1 s s 


(2a) 


Cheng [2] has shown that a desirable BPF is a tuned-loop filter followed 
by a differentiator.. It has a transfer function given by 

2 

F (s) = - T — - g — 2 “ . (2b). 

s + 25 to S + to 
^o o o 


Comparing Fig, 3(b) with 3(a) yields 


F(s) = A 1 F 1 (s) + A 2 F 2 (s) 


(2c) 


and A = 1. Substituting (2a) and (2b) into (2c) 


F(s) = 


A^(s + a) 


A 2 s 


2 2 
s + 25 to S + to 
00 o 


(A. + A„)s 2 + A. (a + 25 to )s 2 + A n (25 to a + 0 J 2 )s + A- ato 2 
1 2 1 00 1 ^o o o 1 o 


3 2 2 

s + 25 to s + to s 
00 o 


(2d) 


Comparison of Equation (2d) and (la) yields 



6 


a 

. 2 
= A,(i)a 

b = 

0 

o 

1 o 

o 

a l 

= A (2£ co a + a) ) 
loo o 

b l = 

2 

to 

o 

a 2 

= A 1 (a + 2 ? q 03 o ) 

b 2 

2? 60 
0 < 

a 3 

“ A 1 + A 2 

b 3 = 

1 


By applying the operator, Equation (lb) yields the following differential 


equation for the M PLL: 


o$e itrc 


2 . 


b^<j> + b^cf) + b + a^sincj) f a^(f>cos(j> + a 2 (<f>coscf) - <j> sin<|)) 


<•<•« c»*f 


••• • ^ *• • •• 

+ a 3 [ C4> “ <t> )cos(j) - 3(f) ^sin<j>] = b^.0^ + b 2 0 1 + b^O^ 


(2f) 


It is understood that tp and 0^. are functions of time. 




7 


III. SOLUTION OF THE M PLL, DIFFERENTIAL EQUATION 

The input signal will be assumed to be a carrier frequency modulated, 
by a single sinusoid. This corresponds to 

0 1 = f3cos(a) t + A) (3a) 

I . m - 

where (3 = Aoj/o) is the peak input phase deviation, w is the modulating 
m ' m 

frequency, ,and A is a phase shift relative to the phase of cj) . It is not 
unreasonable to expect that is a sinusoid of frequency since 0^ and 
©2 are sinusoidal as long as the loop is locked (tracking 0^). That (p is 
sinusoidal has in fact been observed from analog computer simulations of 
the M PLL. Thus it is assumed ; that the loop is. locked and the phase error 
is given by. 

<j) = P sin(w t) (3b) 

m 

where P is the peak phase error. 

There is no, known exact solution to the non-linear 4th order differ- 
ential equation given by Equation (2f). However, under the above assump-- 
tions the steady-state conditions of the solution may be determined by 
the, Principle of Harmonic Balance [3]. Essentially this method consists 
of adjusting the solution parameters so that fundamental components fit 
the non-linear .equation as best as possible. - 

Substitution of Equation (3a) into the right side of Equation (2f) 
yields 

EXR = 0C sin(oo t + A) - 0B cos(rt + A) 
m m . 

= 0[C cosA + B sinA] sin(w t) + 0[C ; sinA - B cosA] cos(U) t) 

m m 

(3c). 
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where 


C = b„m 3 

and 

^ U 2 V 4 

B = b,co - b.os 

2 m 


1 m 3 m 


(3d) 


and EXR means expansion of the right side of (2f). 

Substitution of Equation (3b) into the left side of Equation (2 f ) 
and using the Bessel Function expansions [7] 


cos(P sin(wt)) = J q (P) + ^[^(P) cos(2tot) + J^(P) cos (4tot) + . . . ] 


(3e) 


sin(P sin(wt)) = 2[J^(P) sin(oJt) + J^P) sin(3oit) + ...] 


yields the following expansion for the left side of (2f) (EXL) 

EXL = 2[J.,(P) sin(o) t) + J« (P) sin(3w t) + [a - a„P 2 u) 2 cos 2 (o) t) 

1 m , 3 m o 2 m m 

+ 3a Q P 2 a) 3 sin(u) t) cos(o) t)] + {j (P) + 2[J„(P) cos(2o) t) 
j m m m o z m 

2 

+ J. (P) cos(4o) t) P[a.,w cos(o) t) - a_u) sin(w t) 

4 m ± m , m Z m m 


3 2 3 

+ a»u) (cos(a) t) + P cos (w t))] 
j m m m 


(3f ) 


By using trigonometric identities in (3f) and collecting common terms 
for cos(a) m t) and sin(W Tr> t), Equation (3f) becomes 

,2 


m 


m 


EXL = P[-b 9 w 2 - a w 2 (J (P) + J 9 (P) (J,(P) + J (P)) + (3J (P) 

Z m 3 m o Z Z 13 4 o 


+ 4 J 9 (P ) + J,(P))} +-a 1 w B {J (P) + J 9 (P)>] cos(w t) 
Z 4 1 m o Z m 

„ „ 2a 


+ P[w (b 9 w - b,} + a 9 u) {-J (P) + (• 
~ i z ,-m o 


m 3 m 


a 9 Pco 
2 m 


2 " J 1 (P) + J 2 (P) 


- PJ.(P)}] sin(w t) + H.0.T 
3 m 


(3g) 
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S I n 2 uu m t J cos 2 w,nt , 
2 2 

where H.0.T. indicates the higher order terms - gin- foo tfe z ens 
sin -3u, in t ' m 

s in to- fe-Y, etc. 
m 


Equating the coefficients of 
(Sh) and (3c) gives 


sin(u) t) and cos (to t) 
m m 


in Equations 


g[C cos (A) + B sin(A)] = K-i 
f3[G sin(A) - B cos (A) ] = 


(3h) 


where 



By using the trigonometric identity 


sin 2 (A) + cos 2 (A) = 1 , 


it is easy to show that (3h) yields 


3 


2 


2 2 
K l + K 2 

C 2 + B 2 


(3 j ) 


and 
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g(K C - K_B) 

cos (A) = - r ±—- ~ — (3k) 

C + B ; 

Equation (3j) is the desired solution since it gives the : relation- 
ship between the ; peak phase error (P) and . the peak input signal pha.se 
deviation ((B) in .terms of the M PLL parameters and the modulating .fre- 
quency . 
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IV. THE M PLL PHASE ERROR CHARACTERISTICS 


A typical frequency response for the M PLL considered in this paper 


is shown in Figure 4(a). This response is for operation in the linear 
region of the M PLL - where sin (f) - <p. It has been shown [2] that: 


r = r ■ + — : 
^cL2 S 2o) 


A 


1 


2? cL 1 “nl 


for low-pass 
filter 


BW 


3db 


2to 




n2 


0) 

o 


L for band- 
it pass filter 


A 




BW 


3db 


= 

o o 


+ A, 


(4a) 



? cL is the closed-loop damping ratio, is the natural frequency, 

BW_ is the 3db bandwidth of the filter, and G is the peak ampli- 
Jdb 

tude of the BPF (occurs at co = W „) relative to the d.c. gain of 

m n2 o 

the LPF. 

A digital computer program was developed to determine values for 3 
for various values of P and M PLL parameters using Equation (3j). A 
series expansion [7] was used to determine .the Bessel Function to 0.01% 
accuracy. 

The results for peak phase error of 1.5 radians . (P = 1.5) using the 
same M ; PLL parameters as in Figure 4(a) is shown in Figure 4(b). Note 
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the rapid decrease in $ as u m increases from to the 3db frequency. 

For values of a) removed from the influence of either the LPF or BPF, 

g = P = 1.5. This is to ..be expected since the M PLL does not track 

modulation outside of the filter response and the phase error equals the 

input phase. Also plotted are values for A^ = 0 (the LPF removed) and,. 

for - 0 (the BPF removed). As can.be seen the BPF has little effect 

on the values for frequencies ; near -the LPF. Likewise the LPF has little 

effect on the values for frequencies near the BPF. The influence of the 

LPF extends to : approximately 10 times the ; filter natural frequency. Thus 
to 0 ~ 

if > 100 ) , the two filters may be assumed not to interact and the 
~TT2 nl 

appropriate equation (either with, A^ = 0 or A^ = 0) may be used. The above 
results, were also observed for values of ,P in the range from 0 to 3 radians. 


By. setting A.. = 0 Equation (3 j ) may be simplified to yield results 

toc-B W 3db 

valid for 0 )^ in the BPF region (provided > 10^1^* Getting A^ = 0 in 
Equations (2e) and substituting these values into Equations (3d), (31) and 


then -into Equation (3j) gives 



(4b) 


(4c) 


and 
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K(P) = [J q (P) + J 2 (P) - {J x (P) + 3 3 (P)} + (3 J q (P) • 

+ 4J 2 (P) + J (P)}]- (4d) 

Equation (4d) may be written in terms of J^(P) by using the recur- 
rence, formulas for the Bessel Functions [7], 

(4e) 

From Equations (4a) it is. easily shown, that . 




o 


(4f ) 



Figure 5 shows typical plot of (3 vsi P as obtained from Equation 

(4b). The domain of P is restricted to 0 - r radians since P > 7T causes 

sin P to be negative which in turn causes positive feedback in the M PLL : 

and resulting loss of lock., Here we are interested in the characteristics 

during and before loss of lock. It is interesting to observe the behavior 

of P as (3 is increased for Dw = 0 (u> = co ) . For small' f3, P increases 

mo 

linearly; however, a point is reached (at dB/dP = 0) for which a further 
increase in $ does not -satisfy the equation for any.P near this point. 
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At this point P "Jumps" to some value greater : than tt and results in loss 
of lock. . This is the phenomenon observed by Carden, Thompson .and Cheng 
[2]. It should also be noted that for Eh) far enough removed from 0, the 
"Jump" phenomenon does not occur - the M PLL loses lock without a Jump, in 
P. This can, be more readily seen by referring to Figure 6. Data for 
this figure was obtained from Figure 5 and plotted for constant values 
of 8. For j Dto | > 0.08, and increase in $ (holding Da) constant) is accom- 
panied by an increase in P (in a non-linear manner) until the loop breaks 
lock. For |Dw| < 0.08, an increase in 8 initially results in an increase 
in P until a point indicated by the dashed line is reached. At this point 
a further increase in 8 results in a "Jump" in P and , loss of lock. The 
point at which the "Jump" occurs can.be determined analytically by finding 
dg/dP from Equation (4b) and, setting this value equal to zero. The results 
of this procedure yields 


.... . d V P > 

2[4t o “TT - 

o 


03 


X <°“ 2 + 4 ? o > 


(4h) 


The locus of points obtained from this equation is also plotted in 
Figure 6 and, corresponds to the points at which the slope of the con- 
stant 8 lines is infinite. It should be noted that for a different 4 

o 

and G, Figure 5 and Figure 6 will have the same form but different 
numerical values. The figures of course can be obtained by plotting 
Equation (4g). 
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VI. COMPARISON WITH ANALOG COMPUTER RESULTS 

An Analog Computer Simulation of the M PLL was set up on an EAI 680 

Analog Computer. 8 was increased slowly and the resulting phase error 

was observed. Points were recorded when P equaled 1.5 and when the M PLL 

lost lock. This was done for various 03 near u) and for BPF parameters 

m o 

0) = 10, BW =1, G = .85; 03 = 10, BW = 2,,Q = .85; 03 = 10, BW = 4, 

o o o 

G = .85; and o) = 5, BW = 2.0, G = .95. The LPF BW OJ , = .1 and . £ = .707 

O JQD 

for all cases.. The phenomenon of "Jumps" was observed as the M PLL lost 

lock for 03 near 03 . 

m o 

The results of the analog computer simulation are plotted in Figures 

7(a) through 7 (d) along with the values predicted by Equation (4g). The 

values for loop unlock using the equation were obtained from a plot similar 

to that of Figure 5. For D03 near 03 q (where "Jumps" occur) 8 was the value 

at which d(3/dP = 0. For the other 003 (where "Jumps" do not occur) 8 was . 

the value which, corresponds to a P. of 2.5 radians. This was found to be 

the phase error for which the M PLL lost lock for to where there are no 

m 

jumps. As can be seen, the results agree to within about 10% for. all 
cases,.; The analog computer readings wer.e never greater than the values 
predicted by Equation, (4g). This may be partially due to noise associated 
with the analog computer components.. The assumption of neglecting the 
higher harmonic components, in the solution to the differential equation 
is therefore not very restrictive (at least for the cases considered here). 

Equation (4g) indicates that the loop unlock characteristics are a 
function of ( 0 ^/( 0 ^ (actually Doo) if G (peak gain of BPF) and (open -loop 
damping ratio associated with BPF) remain the same. Normally G is made 
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as close to 1 as possible so that the output voltage in the BPF region is 

well above the noise level. Thus a plot (with w normalized by u ) for 

m n 

u) = 5, BW =1 gives the same results as a plot for those values doubled 
o 

(u).. =10, BW = 2) provided the same G.and i is maintained. This applies 
to the results indicated in Figure 6 as well, since those results were 
derived from Equation (4g). Note that the Q of the filter is the same 
in both cases indicating that the M PLL characteristics are related to 
the. energy stored in the filter. 
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VII.- - CONCLUSION 

The derived equations predict that a "Jump" Phenomenon exists for 

0 ) near the center of the BEF region. The "Jump." in phase is larger for 
m 

w closer to the center of -frequency. Numerical values agree to within 

m - ' 

10% of values obtained from analog computer simulations; 

In order to avoid the "Jumps" in phase (and resulting loss of. lock) 
an operational and design constraint - is that 8 be kept well below the 
values at which this phenomenon occurs. It was observed from analog com- 
puter simulations that negligible distortion of. the . signal occurs for 
P < 1.5. Thus a desirable region of operation for the M PLL is that region 
below the P = 1.5 line in Figures 7(a), 7(b), 7(c) and 7(d). From these 
figures it is also evident .that the subcarrier must be accurately centered 
in the BPF region since the maximum allowable 8 decreases rapidly as is. 
increased or decreased toward the 3db frequencies. 

If the center frequency of the BPf and the 3db frequency of the LPF 
are at least one decade apart the filters do not interact. Equation (4g) 
and its associated plots may then be conveniently used in the design of a 
M PLL. Assuming that and G remain unchanged from design to design. 
Equation (4g) may be used to generate a single plot of 8 vs. Do) to be 
used in the design process. BPF's with the same Q (u^/BW^^) have identical 
characteristics when ..plotted as a function of Deo. 
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Fig. 1. Typical Apollo unified S-band spectrum (baseband) with 
frequency responses of a PLL and M PLL used for its 
demodulation. ^ 






















/3 which causes P« 1.5 
(from equation [3] with 
correct volues for A, a 







Equation 



Fig. 6. Typical plot of Deo vs. P. £ = 0.015, to = 10.0, 









Analog computer for unlock 






Analog computer for unlock 
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